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11th ORDER RESONANCE TERMS IN 
THE GEOPOTENTIAL FROM THE ORBIT 
OF VANGUARD 3 

C. A. Wagner 
ABSTRACT 

The orbit of Vanguard 3 (1959-7A) is strongly resonant with 11th order and 
odd degree terms in the geopotential. It affords an excellent opportunity to 
determine a significant linear constraint between these terms. Tracking data 
on this satellite (in the form of mean Kepler elements) are analyzed over a 
3 1/2 year period in the early 1960’s, which ends with the orbit having just 
passed through perfect commensurability. The eccentricity and inclination 
show the deep resonance variations (up to 2 x 10” 4 in e and .02° in i) with great 
clarity. Previous and current geopotential solutions fail to explain these per- 
turbations. The analysis determined the following constraint for the deep res- 
onant terms (in fully normalized harmonics): 

107 (C S)^ n = (8.6 ±0.2, 7.8 ±0.2) = (C, S) n<n - 6.7(C, S) 13 n 

+ 16.2(C, S) 15 u - 22.9(C, S) 17 u + 19.5(C, S) 19>n - 6.5(C, S> 21<1 

- 7.2(C, S) 23 „ + 1 1.3(C, S) 25fU - 5.0(C, S) 27> n - 3.6(C, S) 29fU +- 

If these 11th order coefficients are of the order of 10~ s At 2 , there will be a 
significant contribution to this constraint for terms at least as high as (25, 11). 
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11th ORDER RESONANCE TERMS IN 
THE GEOPOTENTIAL FROM THE ORBIT 
OF VANGUARD 3 


INTRODUCTION 

The subject of orbital resonance and its application to the determination of 
the geopotential has had extensive discussion in the literature (i. e., Cook, 1961; 
Morando, 1962; Blitzer, 1962; Allan, 1965 & 1967; Kaula, 1966; Kozai, 1966; 
Garfinkel, 1966; Gedeon, 1969; Gaposhkin & Lambeck, 1970, and King Hele, 1972a, 
among many, many other authors and papers). In fact a whole symposium was 
held on the subject in January, 1968 at TRW Systems Group, Redondo Beach, 
California. At that time there was a hope that the longitude dependent geopotential 
(to a fairly high degree) could be determined completely from the strong long 
term deviations of commensurate orbits. Indeed, there were a number of pro- 
posals (i. e., Greene, 1967; Gedeon, Private Communication, 1968) then to 
launch one or a series of special satellites to sample a large number of resonant 
orbits for just this purpose. But it was soon recognized that the magnitude of 
such a program was prohibitive especially in view of the many deep and near- 
to-resonant orbits already available for analysis (Wagner & Douglas, 1969). 

Indeed by the early 1970's the gravity harmonics of order 12, 13 and 14 (and 
past degree 20) had been so determined from near resonant, close orbits 
[Gaposchkin and Lambeck, 1970; Lerch et. al. , 1972a] . The deviations of the 
distant communications satellites in deep resonant (commensurate) orbits have 
served as excellent independent proof that the low order and degree terms of 
recent comprehensive solutions (from precise orbit tracking) have considerable 


accuracy [Wagner, 1972a]. 
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The possibilities of exploiting existing resonant orbits seemed exhausted 
by 1971. However, as Gaposchkin and Lambeck [1970, p. 20] noted, there 
does not exist a general theory of orbital resonance, and this should have been 
a warning that these possibilities were far from closed. Indeed, examination 
of the orbital tracking data discarded by the Smithsonian Astrophysical 
Observatory during preparation of their Standard Earth 2 (SAO SE 2) 
[Gaposchkin, 1970; Gaposchkin and Lambeck, 1970] reveals (with a shock) that 
all were orbits suffering from fairly strong or deep resonance effects. These 
were the 9th order resonant orbit of Telstar 1, and the two 11th order deeply 
resonant orbits of Vanguard 2 Rocket and Vanguard 3. It is ironic that the 
actual solution of the long period resonance problems in these orbits was not 
(as elaborately discussed in Gaposchkin and Lambeck, 1970) through selective 
editing of analytic resonant terms. The solution was to eliminate the offensive 
data altogether! 

The other line of evidence that the possibilities are not closed stems from 
the remarkable new investigations by King Hele’s group (i. e., Gooding, 1971a; 
King-Hele, 1972b) of 15th order resonances on close satellite orbits dragged 
(by the atmosphere) slowly through commensurabilities. These new effects 
were originally observed as an incidental mysterious perturbation in the incli- 
nation evolution of the orbit of Ariel 3, studied to reveal the "super-rotation" 
of the upper atmosphere. 

The observations were of a seeming step decline in the inclination. On 
closer examination the evolution was actually a continuous oscillation of in- 
creasing amplitude and period, finally "reflected" at a critical point, not at 
mid cycle, into an oscillation of decreasing amplitude and period. The critical 
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point (of reflection) was found to be at the time of exact commensurability for 
the orbit with 15th order terms in the geopotential [Gooding, 1971b], R. R. 

Allan, 1971 gives the full explanation for this curious unexpected behaviour as 
an example of an orbit being dragged through a resonance. When the drag forces 
are constant and dominate the resonance forces along track during the passage 
through commensurability, there is an exact solution for the inclination change 
in terms of Fresnel integrals. This solution displays the reflection character- 
istics and the building and decaying oscillations observed in the Ariel data. 

Using the observed inclination changes together with the observations of the 
satellite's longitude, Gooding (1971b) derived a linear constraint among all the 
15th order and odd order terms in the geopotential from the orbit of Ariel 3. 

His simple theory was Allan's (1971) formulation of the appropriate inclination 
rates due to the resonant geopotential terms. (Kaula's [1966] formulation for 
the geopotential rates is completely equivalent and will be used in what follows.) 
Because the orbit was nearly circular only one linear constraint could be easily 
derived. But the way was now open for deriving many other such constraints 
from close orbits observed during resonance passage (King-Hele, 1972b). 

i 

In theory every satellite is a candidate for an infinite number of such passage 
events which could be used for geopotentia.1 determination. In practice though, 
the far satellites will decay too slowly to reach these points in a reasonable time, 
while the close ones will pass through them too fast to be of much use. However 
there is still a large number of medium altitude satellites, with perigees between 
about 500 and 600 km whose orbit decay is just right for observing at least one 
such strong event within 5 to 10 years from launch. And for the more eccentric 
orbits with these perigees, there is the possibility of observing two or more series 
of resonances in each passage. 
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The orbit of Vanguard 3 and Vanguard 2 rocket [n (the 2 body mean motion) 

= 1 1 r evo lut ion s /day , e (orbit eccentricity) = .19, perigee height = 520 km] is of 
this multiple resonance type. These Vanguard orbits offer the only opportunity to 
see significant 11th order geopotential effects similar to the unique opportunities 
afforded the 15th order dragged resonance orbits. But this opportunity has not 
been properly exploited. Indeed the Smithsonian has abandoned them altogether. 
Thus it can be fairly said that 11th order terms, until now, remain among the 
most poorly determined in the geopotential (see Discussion). Unfortunately, 
Vanguard 3 has only been well observed through the first strong resonance (in 
1963-64) which was with odd degree terms proportional to e 2 and thus relatively 
weak. In 1966-72 however, its orbit passed through 3 potentially stronger 
resonance series with even and odd degree terms (still 11th order) proportional 
to e and 1. And in 1973-74 it should be passing through another strong resonance 
(similar to the 1963-64 commensurability) . 

But all of this history is fine hindsight. Even the first resonance pass was 
unknown to me when I examined a set of 182 mean Kepler elements for Vanguard 3, 
calculated by T. Heuring (private communication, 1972) from the Smithsonian's 
precision reduced Baker-Nunn camera observations (Table 1). I knew Gaposchkin 
and Lambeck’s (1970) difficulties with this orbit and that it was strongly affected 
by drag (especially during the solar cycle high of the very early 1960's). My main 
concern was whether the orbit program I was using [see: Wagner and Douglas, 
1970] could handle the strong and variable drag on a 3 1/2 year trajectory. This 
program, ROAD (Rapid Orbit Analysis and Determination), integrates orbits semi- 
numerically using model atmospheres due (originally) to Jacchia, 1965 and 1971. In 
the integrator only terms not containingthe satellite's mean anomaly are retained. 
In this sense the program generates mean elements under a wide variety of forces. 
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But it also can employ empirically derived secular rates and accelerations of 
the elements (to the 5th power). These were found necessary to keep the along- 
track deviations of Vanguard 3 to less than 2° over the long arc when integrating 
with either model atmospheres. (At this level, the longitude of the computed 
satellite is sufficiently accurate to obtain realistic resonant geopotential 
terms.) After a considerable amount of orbital experimentation in ROAD, 
a 3 1/2 year nonre sonant trajectory (with drag and radiation effects) was computed 
for Vanguard 3, with less than 2° along-track variation from the observed 
elements. This trajectory used the Smithsonian Standard Earth 2 zonal geo- 
potential coefficients. My original hope was to "see" the resonance in the 
drag affected mean anomaly observations, but the errors in the model atmos- 
phere were too great to allow this. However the comparison of the inclination 
of this computed trajectory (only slightly affected by drag) with the "observed" 
mean elements showed a clear but unexplained oscillation of increasing am- 
plitude, breaking off near the middle of a cycle at the end of the data (Figure 1). 

This oscillation-breakoff was, of course, reminiscent of the dragged-resonance 
phenomenon so well explained by Allan (1971)for near circular, uniformly dragged 
orbits. Could these residuals be due to the same phenomenon on an eccentric 
orbit subject to a more complex drag-resonance regime? The answer was 
supplied by simulating a Vanguard 3 resonant trajectory with all relevant terms 
for the geopotential harmonics (11, 11) and (12, 11) (from the Standard Earth 2) 
and comparing it with the non resonant trajectory computed previously. The 
form of the inclination residuals produced (Figure 2) emphatically showed that 
the resonant terms would account for the phenomenon. It remained to identify 
and solve for the harmonic constraint(s) that would reproduce the observed 
inclination variations. 
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ANALYSIS 


Kaula [1966, p. 49] identifies geopotential resonance as occuring when the 
orbit longitude ’A is stationary with respect to a particular gravitational harmonic 
term (t , m, p, q): thus 


= 0 = (l - 2p) a + (l - 2p + 1) M + m(0 - 9), 


where and 0. are the orbit's argument of perigee, mean anomaly and right 


ascension of the ascending node and 6 is the rotation rate of the earth. For 
Vanguard 3, the mean motion (m) is close to 11 revolutions /day, which, since 


andfi are comparatively small, establishes m = 11 as the lowest order resonant 
term (under the specification -i- 2p + q = 1). Rewriting (1) for these terms, 
leaves only the q index unspecified: 


'/'res. = 0 = + M) - + ll(n - 0), 

from which; 

q(res. ) = [11(H - 9) + (u + M)] /w • (2) 

From Figure 1, commensurability occurs in the neighborhood of 38400 MJD at 
which time the observations show M = 3992.4 degrees/day, co = 4.9 degrees/day 
and n = 3.3 degrees/day. The solution of (2) (using d- 360*985 degrees/day) 
is very close to q (res.) = -2. Since - 2p + q = 1, then i = 3 + 2p and ■i is odd. 
Thus the resonant series for this commensurability passage is: (^,m, p, q) = 

(' £ , 11, ('&- 3)/2, -2),"C = 11, 13, 15, . . . The resonant longitude i/; for this 
series is simply 

'/'res. = 3w + M + ll(fi- 6) (3) 

The observed variation of \b (= 0 . ) is shown at the bottom of Figure 1 

r res. r ll,~2 

and the stationary value is indeed near the break in the last cycle (^3 8360 MJD). 
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And since it does become stationary, I chose to base the constraint on the 
determination of the inclination rate (I) which is always well behaved during the 
passage. This rate (due the geopotential) is given from the Lagrange planetary 
equations [Kaula, 1966, p. 29] as: 

^ _ cos I BT 1 BT 

na 2 (l - e 2 ) 1/2 sin I ^ na 2 (l - e 2 ) 1/2 sin I ^ 



'hmpq ~ ~ 2 P) w - 2p + q) M + - 0). 

In the above expressions,^ is the earth's Gaussian gravity constant, a e is its 

mean equatorial radius and a is the orbit's semimajor axis. The f functions are 
sinusoidal and the G functions are generally monotonic of order e l q l [see Kaula, 
1966, p. 37]. The C^ m and S^ m are the usual gravitational harmonic coefficients. 
For the resonant series defined above, l-mis even, and \p ( = ip res ) is the same 
for each term in the series. The rate due to the combination of these terms is 
then found to be simply: 
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where 




l 

Ma F. 


I, 11, (l~3)/2 (1-3)2. ■ 

^+3 


[11 - 3 cos I] . 


( 5 ) 


na 3 sin 1(1 - e 2 ) 1/2 
Clearly, only the sum of the sine and cosine terms (in brackets) in equation (4) 

are what can actually be determined from observations of the variations of the 


inclination with (the slowly changing) orbit longitude. 

Calling these terms the lumped coefficients of the 11th order C n and S n , 
the two constraints which can be found from the I data are: 


<c n . v = >n k). ' 

Evaluating the from (5) for the average Vanguard 3 orbit in the resonance 
pass (a = 1.332 earth radii, e = 119, I = 33.3°; the fi/n is irrelevant since it 
does not depend on -i), and normalizing with respect to the (11, 11) term gives 




Qvnlim'f ntnv w-i i n n K 1 a nAM vm-inf n o. 

wvkvx iiiiiiu-tjxt/ rollon aim- ao. 


(C, S) n = (C, S)„ „ - 6.7(C, S) 13 ,, + 16. 2(C, S),, „ 

- 22 « C ' * 19 5(C. S) B u - 6. 5(C, S) 2 , „ - 7.2(C. S),, „ 

t 11.3 (C,S) m n -5.0(C,S) n n -3.6(C,S) M „ .... , (7) 


where the (C, S) u of (6) and the (C u , S n ) of (5) differ only by an irrelevant 
constant factor. 

At this point Gooding, (1971b) solved directly for the lumped co- 
efficients from sets of differenced I values (with observed values of dj ) 

' “res. 

for Ariel 3, using a least squares fit to Equation (4). I chose a more hazardous, 
though potentially more satisfying method to find these coefficients. The method 
was to solve for a set of (15, 11) coefficients from the same program (ROAD) 
which revealed the resonance effects (Figure 1). The hazard taken was in the 
attempt to use all the data in this determination. As I mentioned earlier, the 


8 



drag, model does have deficiencies^ which,>.eyen ■solving for; empirical secular .... 
terms as j high;as the. fifth power --in the mean anomaly.vjstill leaves thercomputed. 
or Jbit' longitude with -peak errors .of as much, as 1 .1/2 degrees. But this level t 
appears) adequate for- a decent discriminatipn of- these- coefficients. ./The -.advantage 
;pf the method ,- of course , is more complete use of the tracking data. In particular, 

dragged-respnance yariatipns .as the., 3 ? 
inclination (Figure 3) v;5 This variatipn'addsito^.the-stren^h of.the constraint 
determination set up for .the inclination data [Equation, (7)] ..But it.-cannot, be 4.. 

used to, determine an independent constraint for these terms. In fact the ec- t , 

■ u : . : mv s.n:‘jv b or,.; m nouo.o&s 6 * vzs'ia'i&'zu o- iv* ?>&*>'•' J>:L 

centricity and inclination variations from these resonant terms are in phase 

a . /on a e: 1 j e.r» *« *o* roo/Bor i; •. oe-'o x . . 

and differ only r by f a constant orbit-factor; namely [from, Kaula, 1966, p. 40] 

r b-e^sVn’l ^ 

flow ca ; W, cps-j,- i,l] ; . fofcw4 ci ^ 

:; This-ratiois-0. ; 663 for- Vanguard 3 * C o mpar is on of the variations 'in Figures 
(1) and' (3) gives the same ratio from^tlie' actual data/ ' ;; jfc*f§iUA w*tH <•• ' • 

The solution (in ROAD) for (15, 11) used the'-S^afidard'EaKh^ w1fch l &ilits iyj,i 
11th order terms [through (16, 11)] having significant long term effect onthe,,^ r . 
Vanguard 3 orbit. Besides the deep resonant series specified by.q= -2 

• KtiiO *i O. ‘i .UP v_ i. *0* J >3 -.OTti 00 a; X) '/.t/i ''.0 0/ "'T-iV 1 . £.•. 

[ adjusted by the solution for (15, 11).]. , the side band series .with q = -1 and -3 

~:.:nw:\ /. sm y.u: re.:. > v.i v.-.«?or> ;.c psai'i-aaeu:) a-,€?oouq uorr.'srroo 

was included; namely the terms (12, 11, 5, -1), (12, 11, 4, -3), (14, 11, 6,, -1), 

tAou.uiOB .u/ro s. -yxJZ! 

(14, 11, 5, -3), (16, 11, 7, -1) and (16, 11, 6, -3). These terms produced effects 

: ' K> .(10.0 < k \ y /0-G ± srt.o; - ( tl . w ? p)'oi 

with a maximum period of 165 days. The weights for all the "observed" elements 

in the* (hffer&ftiaT fcorrfefc't'idn 

cont^i^^ui^'dsplv-ieli^eMe'ct^'dUe'^Wdst likely) tbfdragiwas! ddwhweighted^sb'''as)[^ 
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to influence the (15, 11) solution as little as possible. A similar downweighting 
for the co and H data was found to be necessary because of unresolved systematic 
errors in these arguments. The cause is probably drag again, but not directly. 

The computed semimajor axis shows (even in the final solution) a quadratic 
deviation from the observations of up to 100 m. over the 3 1/2 year data span. This 
is enough to explain the "cubic" fi and oo residuals by way of errors in the computed 
secular rates of these arguments due to oblateness. But even with these model 
deficiencies, the least squares fit of the (15, 11) coefficients to the element 
data produced a dramatic reduction in the residuals for the inclination and 
eccentricity. The overall results of the differential correction is shown in 
Table 2. Also shown in this table are comparison solutions with the unadjusted 
resonance terms of the Standard Earth 2 and the Goddard Space Flight Center's 
(Goddard Earth Model) Gem 4 Field LLerch et. al., 1972b] as well as with the 
original non resonant field. The Gem 4 included data from the deeply resonant 
Vanguard 3 and Vanguard 2 rocket body orbits. The significance of these test 
results will be discussed shortly. 

RESULTS 

The normalized harmonics found from the straightforward ROAD differential 
correction process (described in more detail by Wagner, (1972b) for a multi- 
arc zonal solution) was: 

107 ( C is,u- s i 5> n) = (0.32 ±0.01, °.73 ±0.01), (8) 

with a correlation coefficient of -0.51 between these terms. I was pleased to 
find that the correlations of (15, 11) with the empirical terms in the mean anom- 
aly were all less than | 0.1 1 . This poorly determined data was apparently not 
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influencing the solution because of its deliberate downweighting. But the real 
proof of these speculations, indeed of the (too simple?) constraint theory itself, 
must be to derive a result from that theory and then test it against the data. 

Evaluation of (7) with these adjusted (15, 11) terms and the (11, 11) and 
(13, 11) terms in the Standard Earth 2 (included in the trajectory) yields the 
following values for the Vanguard 3 resonance constraint: 

10 7 (C, S) n = (8.6 + 0.2, 7.8 ± 0.2) ( 9 ) 

Equations (7) and (9) provide the means to obtain a single additional pair of 
coefficients for a field. For example, the constraint using the Standard Earth 2 
field requires (17, 11) to be: 

i° 7 ( c , S) n n = (- .457, .149), 

which appears quite reasonable according to the simple rule proposed by 
Kaula, (1966) [(C, S)^ ^ = 10~ 5 /-t 2 ] . However, as can be seen from (7) this 
result is undoubtedly an accident because the effects of higher degree terms 
are far from negligible. Nevertheless this coefficient set is consistent with 
the Vanguard constraint and should fully explain the eccentricity and inclination 
residuals. Indeed the trajectory calculated in ROAD with this (17, 11) set was 
the equivalent of that calculated with the (15, 11) adjusted values, as shown also 
in Table 2. While this result was expected, it was still enormously gratifying 
as an elaborate numerical confirmation of the simple constraint. 

DISCUSSION 

The results of the previous section can best be summarized in a C, S 
diagram for the 11th order constraint (Figure 4) which emphasizes the relative 
superiority of the SAO SE 2, over the later Gem 4 solution. Clearly neither 
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solution comes close to explaining the resonant eccentricity and inclination 
variations (see also Table 2). But why does the SAO SE 2 field, which does not 
include any resonant Vanguard data, provide a better fit than the Gem 4 which 
does? It is not so surprising when it is recalled that the 7 day arcs of the 
Gem 4 solution would see very little of the variation in these long term effects. 
Gaposchkin and Lambeck (1971, p. 21) felt that analytic term selection might 
solve this problem in their approach. (The idea is to use only those terms 
which have periods of the order of the arc length: inherently impossible with 
straight numerical integration as used for the orbits in the Gem’s). It is inter- 
esting that an unpublished SAO 1969 satellite-only solution called B6.1, with 
11th order resonance terms through 16, 11 and using Vanguard 3 data, does 
come closer to the mean element constraint than any other recent SAO or 
Goddard field. The Smithsonian analytic (really semi-analytic) solutions use 
30 day arcs where possible. This is clearly an advantage where long period 
resonance is encountered. But apparantly it is not sufficient where some terms 
are considerably longer than a month, (as in this case). On the other hand, in 
spite of the fair success of the B6.1 solution, it is possible that SAO abandoned 
this data only because the solutions with it were so highly correlated that the 
results were unrealistic. But that result might have been anticipated and a priori 
constraints applied to limit the coefficient solutions. My own feeling is that since 
it is never a matter of just one term with a very long period, but a spectrum of 
terms which cannot all be absorbed properly, the resulting solution will tend to 
be both unrealistic and highly correlated if only the long period terms are ignored. 

Another comparison solution, Gem 2, (Lerch et al., 1972a), is also shown in 
Figure 4. It differs from Gem 4 only in that it lacks electronic -satellite data. 

The impact on the Vanguard constraint is obviously not significant. 
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What is significant is the probable impact of the odd degree terms beyond 
the 15th. Using the Kaula rule to estimate the size. of the neglected coefficients, 
and taking the root sum of squares of, the estimated terms , [in Equation (7)] . 

gives the concentric circles (of uncertainty) in Figure 4. These circles, can be 
regarded in two ways.- On the one hand they say that the recent comprehensive 
geopotential solutions (without . Vanguard 3 data) can not be expected to reproduce 
the Vanguard constraint because of likely contributions beyond those solved, 
for in these models. ; On the other, hand .they show that solutions, for coefficients 
beyond about jthe 25th degree will benefit strongly from this constraint. 

Of course, even for the less c omplete, solutions, the constraint could still . 
be used (as in the previous section) to satisfy a. ;set,p.f, "lumped" coefficients 
tailored for the i. Vanguard, 3 orbit. ^ r. *. , ... ~, y . ? .. , 

, Beyond the ; llth ; , order commensurabilities, associated with this resonance 
pass (in 1963-64) are series. .for m,=,22, (associated with q = . -4) , m. = 3.3. (q = -6), 
and so on. But the lowest degree's for these are higher (22 and 33. (etc.)) and the 
resonant variations are proportional to e 4 and e 6 (etc.). Thus it can be ex- 
pected that their effects will be slight compared to the 11th order series. In 
fact I have made a solution for (22, 22) from this data and the result was an 
insignificant determination with no appreciable change in the residuals. 

However, the effects of the ev,en w degree "side band" resonances of 11th 
order (terms with,q.= -riband „-3), cannot,.be so, easily dismissed. Their frequen- 
cies are di stinct > enpugh. throughout •.the. .pass to be. able to discriminate two sets 
of even order terms.- Unfortunately the data.is heavily weighted on the "high" 
side of the commensurabiUty. , This, catches, .the longer, (and shortening) period, 



but lesser effects of the (previously commensurate) q = -3 terms. From MJD 
37315 to 38602 (the span of the best data), the period of these effects falls from 
165 to 75 days. The period of the potentially stronger q = -1 terms (commen- 
surable later, in 1966) lengthens from only 45 to 82 days during this time span. 
Preliminary determinations of these constraints are promising but so far in- 
conclusive. At any rate, these determinations are essentially uncorrelated with 
and do not significantly disturb the q = -2 constraint. 

It is anticipated that a fairly strong constraint for even degree (11th order) 
terms will be determinable by the new mean element method from the q = -3 
commensurability which was well observed (by precise Baker Nunn tracking) in 
1959-1960. Unfortunately, only field reduced optical and X-band North American 
Air Defense Command Tracking is available to observe the probably much 
stronger q = -1, 0, and +1 commensurabilities in the late 1960's and early 1970's. 
But their greater strength should provide other sets of independent constraints 
for the even and odd degree terms. 

CONCLUSIONS 

A strong constraint has been found for 11th order and odd degree terms in 
the geopotential which have previously been poorly observed. This has been 
accomplished from analysis of mean element "observations" of the deep resonant 
orbit of Vanguard 3 in 1961-1964. The constraint is not satisfied by recent 
comprehensive geopotential solutions, some of which have used (improperly) 
Vanguard 3 tracking data. However the constraint is undoubtedly more strongly 
affected by terms above the 15th degree than those below, the latter being the 
only ones so far determined in these solutions. Because of this, the determined 
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condition equations (probably) will only be useful in achieving accurate geopotential 
results when solutions for terms as high as the 25th degree are attempted. 

Examination of independent 11th order commensurabilities for Vanguard 3 
in 1959-60, 1966-71 and 1973-74 (using mean elements derived from tracking 
data) should provide additional strong constraints for 11th order harmonics. 
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92.7024 
— 55*7253 

29.0994 

4.391-1 

341 .88 79 

-34 6.44-96 

289.6278 

- 2 £2 . 058 4 

203. 7764 

- 155*4157 

1 12.181 e 

74*8-950 

43. 191 4^ 

1 - 7 * 0600 - 

356. I77e 

-330. 1366 

304. 1323 

-2 72*2903 

229.7370 

—161 *5936 

133.7755 

— 90 * 447-7 

59.0449 

28*5724 

7.2560 

-3 42.0 530- 

3 17.0539 

286*3385 

249.6067 X 

207*7694 

155.3100 

- 114.2025 

73.0303 

- 43.6456 

16.1249 X 
356.2871 
329.3339 


Table 1 (Continued) 



TIME (MJD) 

A(e.r.) 

E 

INCL (°) 

OMEGA (°) 

NODE (°) 

MEAN (°) 

383 64.C69400C0 

1 . 33 2242 9C 

• 1890459 

33 .3416 

75.4874 

234.3596 

3C4.9032 

38371 .0C41 CCCO 

1 . 33 222 £ 9 C 

• 1890407 

33.3396 

109.3499 

21 1 . £5 15 

270 .83 1 2 

38378. 1667CCC0 

1 .33220950 

• 188811 5 

33 .3402 

143.6613 X 

166.4539 X 

231 •7401- X— 

38385. C4660CC0 

1 . 3322086 0 

. 1885578 

33.3457 

178 .0 137 

165.3642 

181 .5490 

38392 • 1597CCC0 

1 .3322 1040 

. 1 882899 

33 .3519 

2 12.8660 

141 *9805 

131*8222 

38399.00280000 

1 .3322062 C 

• 1881985 

33.3564 

246.4287 

119.48E4 

92.0516 

384C6.74 7200CO 

1.332178 70 

. 1881 853 

33.3577 

264 .4 192 

94.0267- - 

51.1773- 

38413 .OE420CC0 

1. 3321936C 

. 1 883145 

33.3567 

3 15.3439 

73.2925 

31.3251 

38430 .C 8I90CC0 

1 # 33 2 1 6 7 7 C 

• 1889216 

33.3470 

38.7237 

17.2992 

3 34.4366— -- 

38424 .045800 C 0 

1.3321 6560 

• 1 890397 

33.3457 

58.0989 

4.2640 

2 20.3200 

3844 1 .0 7080000 

1.33219010 

. 1 890891 

33.3446 

92.3946 

341*1614 

267.8939 - 

38448. 184 1 0CC0 

l. 3321 7060 

. 1 889897 

33.3445 

1 27 .1437 

217.7642 

250.4623 X 

38469.06940000 

1. 3321 649C 

.1882337 

33.3509 

229.4007 

249.0827 

113*6304 

38476, 00280000 

1.332 161 1 C 

• 1 881533 

33.3518 

263.4 100 

226.2646 

75.652C 

38483.0278CCC0 

1 .332 171 90 

• 1681714 

33.3535 

297.881 I 

203 * 1 €64 

43 * 560 8 

38490 • 144400C0 

1. 33216440 

• 1 883088 

33.349 l 

332.7771 

179.7879 

17.5849 

38497.08 1 90 C CO 

1. 33216540 

• l 885342 

33.3445 

6 .780 4 

156.9746 

356.3611 

38504. 19860CC0 

1. 3321602C 

.1887957 

33.3393 

4 1 .6026 

133.5663 

330.5612 

385 1 1 .44440000 

1 .3321496C 

. 1 889608 

33.3365 

75.0650 X 

1 1 1 .0468 X 

303* 533-7- X— 

3851 8. C6940CC0 

1.33215430 

. 1 889299 

33.3392 

109 .3747 

87.9399 

271 .9805 

38525. 09170000 

1 .3321591 0 

. 1 887574 

33*3435 

143.6897 

64.6452 

229.3736- 

38522 . 1 1250CC0 

1 • 2321582C 

.1881 720 

33.3483 

178 .0657 X 

41 .7576 

181 • 1557 

38539. 0431 0000 

1 .332 1561 C 

.1873837 X 

33.3668 X 

2 12. 1216 X 

18.9379 X 

132.401 6 

38546. 1 5 560 C C 0 

1. 3321437C 

. 1880374 

33.3557 

246 .8997 

355.5901 

9C.C69C 

38553 .000C0CCQ 

1 .33215520 

. 1880127 

33.3554 

280.4723 

333.0904 — 

S7.6C32 

38560 .55830000 

1 .3321541 0 

.1881506 

33.3521 

322 .4778 / 

304.9497 X 

28.1247 X 

38567. 41390CC0 

1 .33214070 

• 1883386 

33 .3475 

351 .1806 

285.6970 

7*4338 

36574. 35000CCO 

1 . 3321 4 89 C 

. 1886 46 8 

33.3407 

25. 149 4 

26 2 .69 C 7 

241.0547 

38581 .01670000 

1 .33214330 

• 1 888640 

33.3367 

57 .7366 

240 .962 1 

2 16.9138- 

38588.13 190000 

1 . 332 1 433 C 

. 1 889669 

33.3319 

92.5004 

217.5553 

287. 96e i 

38595. C6530000 

1.3321 439 C 

. 1888354 

33.3340 

126.3759 

194.7475 

250*5997 

38602. 06750C00 

1.33215130 

. 1 885981 

33 .3390 

160 .7 19 1 

171.6499 

208.0961 

STANDARD OBSERVATION 







SIGMAS 

C.OOCOIOCC 

0.0000200 

0.000500 

0 .020000 

0.020000 

1.000000 


*These elements were originally determined {by T. Hearing) as osculating values from a precision tracking program (GEODYN) using the precisely reduced 
Baker-Nunn optical observations of the Smithsonian Astrophysical Observatory’s World Wide Network. The mean elements are these osculating values 
less the short period terms (due to oblateness) due to D. Brouwer, 1959. The aic lengths to which these original elements apply are generally 7 days. 
After processing the 3-1/2 year arc of mean elements it was found necessary to discard only 48 or 4.4% of these observations. The rejected values are 
marked with an X. 
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Table 2 

Results of ROAD- Vanguard 3 Orbit Determinations 
(Using 3-1/2 Years of GSFC Mean Elements*) 


Field 

Overall 
Weighted 
RMS Residual 

Weighted Element Residuals 
(rms) 

Used 

e 

(2 x 10" 5 ) 

I 

(0.0005') 

Ci) 

(0.02°) 

n 

(0.02°) 

M 

(l.°) Unit Wts. 

SAO SE2 
Without 
Resonant 
Terms 

3.17 

2.29 

6.79 

0.75 

1.16 

0.7 

SAO SE2 
With 

Resonant 

Terms 

2.19 

1.44 

4.49 

1.06 

0.89 

0.8 

GEM 4 
With 

Resonant 

Terms 

2.78 

1.83 

5.92 

1.07 

1.03 

0.8 

SAO SE2 
With 

Resonant 
Terms + 
(15, 11) 
Adjusted 
From data 

1.13 

0.76 

1.59 

1.13 

0.73 

0.7 

SAO SE2 
With Res. 
Terms + 
(17, 11) 
Solved 
from 

Constraint 

1.14 

0.73 

1.59 

1.19 

0.75 

0.7 


*AII orbit determinations used all six mean elements as observed data and adjusted the six initial ele- 
ments, a radiation pressure and drag coefficient and 5 mean anomaly "Secular" terms. Included are 
short and long period luni-solar terms to the 4th degree in their disturbing potentials, as well as all 
significant long period zonal variations in the given fields and the effects of precession and nutation 
of the earth's polar axis. J 
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38606 MJD) ARE SAO MEAN 
ELEMENTS FROM FIELD REDUCED 
BAKER-NUNN TRACKING. THIS 
DATA WAS NOT USED 
IN THE COEFFICIENT ANALYSIS. 


RESIDUALS ARE “OBSERVED” 

MEAN ELEMENT VALUES MINUS COMPUTED 
- VALUES FROM A FITTED ROAD TRAJECTORY 
WITHOUT RESONANT TERMS BUT INCLUDING 
_ ADJUSTED DRAG & RADIATION PRESSURE 
COEFFICIENTS, 6 INITIAL EPOCH ELEMENTS 
AND 5 SECULAR RATES (ACCELERATIONS) 

OF THE MEAN ANOMALY. 
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Figure 1. Measurement Residuals In Inclination From A Five Year Nonresonant Vanguard 3 Trajectory 
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RESIDUALS ARE VALUES COMPUTED 
WITH RESONANT EFFECTS OF (1 T, 1 1) 

AND (12, 11) [FROM SAO SE 2] MINUS 
VALUES COMPUTED WltHOUT THESE 
EFFECTS. BOTH COMPUTED TRAJECTORIES 
HAVE THE SAME DRAG AND OBLATENESS 
PERTURBATIONS. THE NON RESONANT TRAJECTORY 
HAS BEEN “FITTED" TO THE RESONANT TRAJECTORY 
(BY THE ROAD LEAST SQUARES PROCESS). 

THE RESONANCE IS AT A DIFFERENT TIME THAN FOR 
THE ACTUAL VANGUARD 3 ORBIT BECAUSE THE 
ENERGY OF THE SIMULATED RESONANT 
TRAJECTORY IS SIGNIFICANTLY DIFFERENT 
FROM THE ACTUAL. 
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Figure 2. Inclination Residuals From Two 5 Year Simulated Vanguard 3 Trajectories 
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Figure 3. Measurement Residuals In Eccentricity From a 5 Year Nonresonant Vanguard 3 Trajectory* 
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CIRCLES INCLUDE 

LIKELY CONTRIBUTIONS TO THE 1 1TH ORDER 
CONSTRAINT FOR ALL TERMS 
EQUAL TO OR GREATER THANi 

© 11TH ORDER VANGUARD 3 
CONSTRAINT [EQUATION (7)] 

(WITH la NOISE VARIATION) 
DETERMINED FROM MEAN ELEMENTS. 

CONSTRAINT CALCULATED FROM 
• SAO SE 2 


CONTAINS 
VANGUARD . 
3 TRACKING 
DATA 


'■SAO B6.1 FIELD (UNPUBLISHED, 1969) 

O GEM 4 
.□GEM 2 

SAO SE 2, GEM 4 AND GEM 2, ARE 
COMBINED SATELLITE-SURFACE GRAVITY 
FIELDS. SAO B6.1 IS A PURELY OPTICAL- 
SATELLITE FIELD. ALL THESE 
FIELDS CONTAIN 11TH ORDER TERMS 
THROUGH 16TH DEGREE. 


Figure 4. Vanguard 3 Constraint From Various Fields 



